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Abstract 

A model is suggested to unify the Einstein GR and Dirac Cosmology. There is one adjusted parameter 
62 in our model. After adjusting the parameter hi in the model by using the supernova data, we have 
calculated the gravitational constant G and the physical quantities of a(t), q(t) and p r (t) / pb(t) by using the 
present day quantities as the initial conditions and found that the equation of state parameter we equals 
to -0.83 , the ratio of the density of the addition creation fi\ = 0.8 and the ratio of the density of the 
matter including multiplication creation, radiation and normal matter J7 m = 0.2 at present. The results are 
self-consistent and in good agreement with present knowledge in cosmology. These results suggest that the 
addition creation and multiplication creation in Dirac cosmology play the role of the dark energy and dark 
matter. 
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I. INTRODUCTION 



According to the Dirac's arguments, the large dimensionless numbers provided by atomic physics 
and astronomy of our universe are connected with each other 1]. These numbers include: (i) the 
ratio of the electric to the gravitational force between an electron and a proton a\ = e /Q mp ra e ~ 
10 39 ; (ii) the age of the universe expressed in terms of the atomic unit 02 = "Ifg ~ 10 39 , where 
H is the Hubble constant; (hi) the mass of the part of the universe which is receding from with 
a velocity v < c/2, expressed in the units of the proton mass, say, 03 ~ 10 78 . Dirac introduced a 
large number hypothesis 

a x 9* a 2 = a 1 / 2 (1) 

Based on this hypothesis, Dirac suggested a model of cosmology with a varying gravitational 
constant G and an increase in the amount of matter in the universe. 

Based on the large number hypothesis, a number of cosmological models with a varying gravi- 
tational constant have been proposed However most of them met many difficulties. Noting 
that gravitational constant G cannot vary in general relativity (GR), the first difficulty is that one 
must explain the contradiction between Einstein GR and Dirac Cosmology. Though much effort, 
for example, Milne two time scale hypothesis [jj, Weyl's geometry Q etc., has been devoted to 
reconcile the requirements of these two theories, it is still an open question to establish a theory 
which can unify the Dirac cosmology and Einstein GR. 

The second difficulty comes from experiments. Almost all experiments at the scales of solar 
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. A possible variation of G has been 



system and galaxies have not found the variation of G 
investigated with no success through geophysical and astronomical observations. From experimen- 
tal results one tends to believe that G is a constant for local system with large scale. 

The third difficulty belongs to the conservation of energy and momentum. Usually we use a 
perfect-fluid energy-momentum tensor to describe the matter of universe and it is conserved in the 
cosmic evolution. The addition creation, the new matter created uniformly in the whole space, 
and the multiplication creation, the new matter created in regions where old matters exist, must 
come from other mechanism as suggested by Dirac. 

Many years ago a possible unified theory of Dirac cosmology and GR to overcome the above 
difficulties was suggested in i|. The basic idea is as follows. Though one would expect a constant 
value of G at the local system such as solar system, binary system and galaxies, it must be stressed 
that the cosmological observations still cannot put strong limits on the time variation of G in the 
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111 ]. Though Einstein 



cosmological scale, especially at the late time of the universe evolution 
GR has been proved to be correct by many experiments such as the excess perihelion precession of 
Mercury, gravitational redshift etc. at the scales of local system, it probably needs to be modified 
in the cosmological scale, especially if we want to use it to explain the recent observational result 
of the acceleration. Using the idea of Dirac cosmology with the variation of G one could provide 
a possible way to modify the Einstein GR. 

As is well known, the Lagrangian density of Einstein GR with a cosmological term is 

where a = 2A/R is a dimensionless parameter. The cosmological constant A can be generally 
explained as the background fluctuation of the cosmological vacuum, and R is the 4-dimensional 
scalar curvature. In ordinary astrophysical problems of local systems, the magnitude of R is about 
8irGp + 4A. Since the vacuum energy density is much smaller compared to the density p in local 
system, a is a very small quantity and can be neglected. But in cosmological problems, because p 
is small and has the same order as that of the vacuum energy density, a cannot be neglected. In 
the limit of vacuum or the matter domination area, p — > 0, R — > 4A, a can attain the magnitude 
1/2. This means that a plays an important role at cosmological scale. However, there is only a first 
order term of a in the equations of Einstein cosmology. Instead of the factor 1 — a, we argue that 
in a perfect cosmological theory, the Lagrangian density of the gravitational field could contain 
higher-order terms of a. We take the Lagrangian density as 

Le = J^^9 R f( a )> (3) 



where 



The Einstein action becomes 



f(a) = 1 - a - b 2 a 2 - ... (4) 



S = J dAx V-gR{ 1 - a- b 2 a 2 - ...) + J d A xyf=gL M - (5) 

Starting from Eq.l(5]). we can establish a Dirac cosmology with varying G and matter creation. 
Obviously, this theory can unify GR and Dirac cosmology because in a local system with large 
scale, a —* and the theory reduces to Einstein GR. Only at cosmological scale, the terms a, a 2 ... 
become important and our theory reduces to Dirac cosmology. 

Recently, many authors have considered the terms R m (m > 1) and/or R n (n < 0) on gravity 
^3]-^^]. But none of them has connected with Dirac cosmology and the large number hypothesis. 



In this paper, after introducing the higher order terms of a, we will establish a unified theory of 
GR and Dirac cosmology and explain the acceleration of the contemporary universe, variation of 
G, cosmological supernova Type la data, and the creation of matter. We find that the addition 
creation of matter looks like the dark energy and the multiplication creation of matter like the 
dark matter. 

The organization of this paper is as follows: we will present our model in Sec. II. In Sec. Ill, by 
using the present cosmological parameters as our initial conditions, we will calculate numerically 
the physical quantities of our universe and study the addition creation and multiplication creation. 
We will prove that our results are in good agreement with present knowledge in cosmology. Finally, 
a summary and conclusion will be given in Sec. IV. 

II. UNIFIED THEORY OF GENERAL RELATIVITY AND DIRAC COSMOLOGY 

Up to the second order of magnitude 0(q 2 ), the action of our model is 

S = J d A Xy/^gR(l - a- b 2 a 2 ) + J d A x^/^g~L M , (6) 

where b 2 is a expansion parameter and we will adjust it from supernova data later. The field 
equation for the metric is then 

(1 + b 2 a 2 )R^ - ^R(l -a- b 2 a 2 ) 9fiu + b 2 a 2 R 2 (g^V°V a - V^V,)ir 2 = SirGT™, (7) 

where Tj£ is the energy-momentum tensor of matter and V^(^ = 0, 1, 2, 3) is the covariant deriva- 
tive induced by the metric. When b 2 7^ 0, the constant curvature vacuum solutions require 

-l±^T+W 2 

a = 3^ ' (8) 

We find the constant-curvature vacuum solutions are not Minkovski space, but rather de Sitter 
space (+) or anti-de Sitter space (— ). 
Eq.Q can be rewritten as 

G^ = %txG{T^ + 9^), (9) 
where G^ u is the Einstein tensor, and 

G = G/(l + b 2 a 2 ) (10) 



A 1 Ub 2 A 2 (G-G) 1 . . 1 

^ = ~8^G 9 ^ ~ 8^GV G ^ " 8^ V Vct " V ^ } G- 
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The physical meaning of G and 9^ v is obvious. We see from Eq.Q that G plays the role of 
gravitational constant and 9^ v the creation of energy-momentum in the universe. They depend on 
the second order term of the expansion f(a). If 62 = 0, G = G = const., 9^ u reduces to the 
ordinary vacuum energy cosmological term. Since G, 9^ v depend on R and then are functions of 
time t, this is just the character predicted by Dirac large number hypothesis. Futhermore, using 
Eq.@ and the Bianchi identity, we get 

= -(T# + e^VhiG. (12) 

9^ v is no longer conserved. To understand the physical meaning of Ea. (|12j) . we rewrite it as 
following: 

V v 6$=-T$V v hiG (13) 



V^ 2 J = -VV^lnG, (14) 

where 9^J +9$ = 9a V - It is clear that 9$ is created around the normal matter, while the creation 

of 9$ spreads over the whole universe. In the words of Dirac Cosmology, 9$ is the multiplication 

(2) 

creation having the basic character of dark matter and 9\{, is the addition creation having the basic 
character of dark energy. Both of them depend on the variation of G. In Sec. Ill, the suggestions 
that addition and multiplication creations being dark energy and dark matter respectively will be 
elaborated numerically. 

Employing flat Robertson- Walker metric 

ds 2 = -dt 2 + a 2 (t)dx 2 (k = 0) (15) 

and the perfect-fluid energy momentum tensor 

= (PM + VMWpPv + PMQimu, (16) 

where W 1 is the 4- velocity in fluid-rest frame, pm, Pm are the energy density and pressure respec- 
tively. Taking pu = wpM as the equation of state, w = for matter, w = 1/3 for radiation, we 
find that the field equations become 

3H 2 = 8ttG(pm + pe) = SttGpm + |v(GQ + 3# (17) 



2H + 3H 2 = -8ttG{ Pm + Pe) = -8ttG Pm + qV(G) - 2^(£ + (18) 
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Figure 1: The potential V(G) with b 2 = 1 

where H = a/a, pe = #oo, pq = da/gn (i = 1,2,3) are the energy density and pressure of the 
created matter, respectively. Comparing to the kinetic term (^) 2 , 



(19) 



which corresponds to a potential. Here we only consider the case of positive curvature scalar R. It 
is easy to show that R is always positive in the evolution of the universe with matter and radiation. 
The kinetic term and the potential determine the evolution of G explicitly. V(G) will reach its 
maximum at G m = G/[l + 62/(1 + \/l + 3&2) 2 ] as shown in Fig.JQ. 

We also define p eW = 6$, p e(2) = 9$, p eW = O^/m and p e{2 ) = 0$ /gu (i = 1,2,3) as 
the energy density and pressure of both of the multiplication and addition creation matter. The 
components of /i = of Eq, ()13|) and Eg. 1)14(1 can be written as 



p ii) +3H(p gW +p $ m) 



G 

~ PM G 



(20) 



p Q {2) + 3H(p e(2) + p e( 2) ) = -p e -= (21) 

The other components of Eq. ([T3|) and Eq.p4|l(// = 1, 2, 3) give the identities about pressures p e (i) 
and pg(2) and no further information. Since the multiplication creation #W looks like dark matter, 
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we set 

Pew = Pe(2) = Pe (22) 
Ea.(|17 j) -Ea.(|22 |) and the equation of scalar curvature 

R = 6[- + (-) 2 ] = 6(H + 2H 2 ) (23) 
a a 

given by the flat Robertson- Walker metric form a complete set. In principle, we can solve the set 
of equations to find the behavior of the cosmological evolution. 

III. THE NUMERICAL SOLUTIONS OF THE UNIVERSE WITH MATTER AND RA- 
DIATION 

We now study the behavior of our universe containing matter and radiation. To find the 
cosmological solution, we need the initial conditions of Eq.(|17 |) -Eq.(|22j l in the early universe which, 
unfortunately, is little known now. Hence, we have done our numerical computation by using the 
present cosmological parameters as our initial conditions, which include the Hubble parameters 
constant Hq = 0.7 x WOkm ■ s _1 Mpc _1 , the present deceleration qo = —0.5, the cosmological 
constant A = 1 (as the scale), the density of baryon Qf, = 0.05, the ratio of the baryon-to-matter 
(including baryon and multiplication creation) £lt,/£l m = 0.17 and the ratio of baryon-to-photon 
r\ = 6.1 x 10 -10 [2]]]. We will start from the present epoch and trace back the history of our universe. 

FigEl shows the variation of G in the evolution. The present time is set as to = 0. The three 
curves correspond to 62 = 0.66, 2.5 and 10 respectively. We see that the 'gravitational constant' G 
increases at first when the universe expands and then decreases at the late time. The ratio of the 
density of radiation to baryon floats up quickly and the scale factor a(t) — ► i 1//2 as t approaches 
the age of the universe, which are shown in FigEJand Fig^l respectively. Therefore, the radiation 
dominates the early universe. And the age of the universe given by a(t) = is 1.02Gyr, 1.16Gyr 
and 1.34Gyr for 62 = 0.66,2.5 and 10 respectively. Since there is no mechanism of inflation in 
our model, it is easy to understand that our model can consider the universe back to the epoch 
of radiation domination. FigO shows the evolution of the deceleration factor q. We see that the 
universe decelerates in the early era and gradually stops decelerating and starts accelerating. This 
is in a good agreement with what we understand the universe nowadays. Through the deceleration 
factor q, we notice that the scalar curvature R is always positive. 

Obviously, 62 is the parameter in our model which needs to be adjusted. To determine this 
parameter, let us compare our model with the supernova data. Noting that the uncertainty of 
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Figure 3: The ratio of the density of radiation-to-baryon p r /pb versus time t. 

the present value of the deceleration factor go is still large, we will give a better estimation of 
the parameter 62 and qo to fit the supernova data. FigEl presents the supernova data as a phase 
portrait of the universe ^| 2^], where the universe was decelerating at high redshift and started 
accelerating when it was about two-third of the present size. Our model is favored by this set of 
data : 62 = 2.5 and qo = —0.5. Comparing FigEJand the curves given by Fig. 2 of ref 20], we find 
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Figure 4: The scale factor a(t) versus time t. The three curves from right to left correspond to hi = 0.66, 2.5 
and 10. All of them are asymptotic to t 1 / 2 as t approaches the age of the universe. 

1.0 
0.5 
0.0 
-0.5 

-1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 
time t unit in 10 10 yr 

Figure 5: The evolution of the deceleration q(t). The three curves from right to left correspond to 6 2 = 
0.66,2.5 and 10. 

that the fit of our curve is better than that of the standard model. 

Fig[7| shows the evolution of the pressure of the density of the addition creation matter with 
&2 = 2.5. The pressure is highly positive in the beginning of the universe, falls quickly blow zero and 
has a negative value in the recent and present era. This negative pressure looks like the character 
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Figure 7: The evolution of the pressure of the addition creation p 6 {?) at 62 = 2.5. 

of the dark energy. If we take p d {2) = wop e (2) as the equation of state of addition creation, we show 
the parameter of the pressure-to-density of the addition creation wg versus time t in Fig|H]and find 
wg(t = 0) = —0.83 which is much smaller than —1/3 and is consistent with the requirement of 
dark energy to explain the accelerating universe. To compare our results with other models, we 
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Figure 8: The evolution of the parameter of the pressure-to-density of the addition creation wg versus time 
t at b 2 = 2.5. 



have also calculated the ratio of the density of addition creation and find 



Po<?) 



= n A = o.8 



(24) 



Pg(l) + Pq{2) + p r + Pb 

and that of the density of matter (multiplication creation and normal matter) in our model reads 

PqW + p r + Pb 



n M = o.2. 



(25) 



Pg(l) + pg(2) + p r + pb 

These results have the same magnitude as that dark energy and dark matter respectively. Hence, 
we suggest that the dark energy comes from the addition creation and the dark matter from the 
multiplication creation. 



IV. SUMMARY AND DISCUSSION 



We have suggested a model to unify the Einstein GR and Dirac Cosmology. In local system, our 
theory reduces to GR, but in the cosmological scale, our theory refers to the Dirac cosmology. The 
variation of the gravitational constant comes from the scalar curvature. In the local system, the 
variation of the gravitational constant is negligible. This result is in good agreement with present 
experiments. But in the cosmological scale, the change of G is remarkable. The acceleration of the 
present universe is a proof of the decrease of the gravitational constant G, because the decrease of 
G corresponds to an effective repulsion. 



11 



After introducing the term —62a 2 in our theory and fit the adjust parameter 62 by supernova 
data, we have calculated the physical quantities of G(t), a(t), q(t), p r (t)/pb(t) by using the data of 
present epoch as the initial conditions. We have found that the results are self-consistent and in 
good agreement with present knowledge of cosmology. 

According to Dirac large number hypothesis, matter will be created in the universe. We have 
calculated the matter comes from addition creation and that from multiplication. An interesting 
picture in our theory is that the addition creation, which spreads over the universe, looks like the 
dark energy and the multiplication creation, which clusters around the normal matter, like the 
dark matter. We have found that the pressure of addition creation has a big positive value initially 
and fall down quickly to become a negative value. It has a negative pressure region corresponding 
> t > —0.74 x 10 10 yr. The equation of state parameter wg = —0.83 at present. This value is in 
agreement with present dark energy models. We have also calculated the ratio of the density of 
the addition creation and found S7a = 0.8. The same parameter but for multiplication creation, 
radiation and normal matter has also been computed, which reads {l m = 0.2. Both of them have 
the same magnitude of the observational value of the dark energy and matter. This result suggests 
that the dark energy and dark matter are just the addition creation and the multiplication creation 
in Dirac cosmology. 

Finally, we would like to emphasize that this model cannot be extended to the big bang epoch. 
We have not add the terms with R m (m > 1) (or oT n (n > 1) ) in the Lagrangian density. Obviously, 
these terms are very important to the very early universe, especially in the inflationary epoch. Our 
model can only be used in the time evolution regions starting from the radiation dominated epoch 
to the present time. 
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